from partie3 import *


#=======resolution du systeme et tracage de la solution========#

I=array([1,0,0])
#y = resoudre(I, 0.1 , 0.1 , 1. , 2.)
#courbe = mp.plot(y[0],y[1])
#mp.legend(courbe ,("X","Y","Z"))
#xlabel("t")
#ylabel("y(t)")
#mp.show()
#======calcul du temps du debut de la phase periodique========#

#y = meth_epsilon(definir_sys(0.12),I,1. ,2.,0.1,"rk4","non")
#print partie_periodique(definir_sys(0.12),0.,y)


b=0.1
eps=0.1
debut=1.
end=1.



#=======Tracage du diagramme de bifurcation======#
B=linspace(0.1,0.2,20)
#bifurcation (I,B,0.1,1.,2.)
